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Abstract 


Land subsidence induced by soil consolidation is one of the natural hazards that occur gradually. This 
phenomenon has reached its critical state in most regions of Iran. Factors affecting subsidence are 
groundwater level depletion, land cover, soil type, elevation, slops gradient, rock bed depth, etc. The 
causes of subsidence should be investigated so that decisions could be based on the real characteristics 
of the region. Decrease in groundwater level is one of the most important factors that influences 
subsidence but due to the complexity of the relationship between subsidence and other factors, a direct 
linear relationship between groundwater level and subsidence cannot be considered. The study aimed to 
investigate the relationship between subsidence and groundwater changes through relationship between 
the two parameters in the period 2014 to 2018 in Isfahan. In this regard, the time series of radar 
interferometry and the time series of water levels of piezometric wells in the Isfahan were studied. The 
results show that in 12% of the wells, no correlation exists, in 9%, poor positive ignorable correlation, in 
12%, fair positive correlation, in 48%, significant positive correlation, in 6%, poor negative ignorable 
correlation, in 4% fair negative correlation and in 9% there is a strong and negative significant correlation 
between the subsidence rate and the groundwater level. Correlation values show the complexity of the 
relationship between subsidence and water level depletion. Therefore, the relationship between 
subsidence and its causative factors requires more detailed studies and comprehensive models. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و چهارم؛ زستان ۱۶۰۱ صص ۱۷۳-۱۹۱ 
مقاله پژوهشی 


ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران يا بحران مدیریت) 


A e Me ۰ E هر‎ 7 5 ae 5 ae 0 >. “fe 
تداخل‌سنجی راداری (منطقه مطالعاتی: اصفهان)‎ 

زهرا آزرم- دانشجوی کارشناسی ارشد سنجش‌ازدور گروه نقشه‌برداری, دانشکده حمل‌ونقل, دانشگاه اصفهان 

® حمید مهرابی '- استادیار گروه مهندسی نقشه‌برداری, دانشکده عمران و حمل‌ونقل. دانشگاه اصفهان. 


سعید نادی- استادیار سابق گروه مهندسی نقشه‌برداری, دانشکده عمران و حمل‌ونقل. دانشگاه اصفهان. محقق دانشکده مهندسی عمران و 
محیط‌زیست. دانشگاه كارلتونء اوتاواء کانادا 


تاريخ دریافت: ۱۶۰۱/۱/۱۵ تاريخ بازنگری: ۱۶۰۱/۳/۱ تاريخ تصویب: ۱۶۰۱/۳/۱۹ 


حکیده 


فرونة ت زمین یکی از مخاطرات طبیعی است که به‌صورت تدریجی اتفاق می‌افتد. اين پدیده در اکثر 
مناطق ایران به وضعیت بحرانی رسیده است. عوامل تأثیرگذار بر فرونة ت شامل تغييرات سطح آب 


رخداد آن بررسی شود تا با توجه به ویژگی‌های هر منطقه تصميم كيرىها مبتنی بر واقعیت باشد. افت 
سطح آب زیرزمینی یکی از عوامل مهم موثر بر فرونشست است ولی به دلیل پیچیدگی ارتباط فرونشست 
با عوامل تأثیرگذار دیگر. نمی‌توان Sy pow‏ عام يك رابطه خطى مستقیم بين تغییرات سطح آب زیرزمینی 


و فرونشست در نظر گرفت. اين تحقیق با هدف بررسی رابطه بين فرونشست و تغییرات سطح آب 
زیرزمینی به بررسی همبستگی بين دو پارامتر مذکور در بازه زمانی ۱۳۹۳ تا ۱۳۹۷ در منطقه اصفهان 
پرداخته است. در اين راستا سری زمانی تداخل سنجی راداری و سری زمانی سطح آب چاه‌های پیزومتری 
در منطقه اصفهان موردبررسی و تحلیل قرار گرفت. نتایج نشان می‌دهد. در ۸۱۲ از چاه‌ها. عدم همبستگی» 


Email: h.mehrabi @eng.ui.ac.ir ۰۳۱۳۷۹۳۵۲۹۷ نویسنده مسئول:‎ ١ 
نحوه ارجاع به اين مقاله:‎ 
بررسى ميزان همبستگی بين فرونشست و تغييرات سطح آب زيرزمينى با استفاده‎ VEN) آزرم» زهرا؛ مهرابی» حميد؛ نادی» سعيد؛‎ 


از آناليز سرى زمانى تداخل‌سنجی رادارى (منطقه مطالعاتى اصفهان). جغرافيا و مخاطرات محيطى. .)٤(۱۱‏ صص ۱۷۳-۱۹۱ 
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۱۷ جغرافیا و مخاطرات محیطی byl‏ چهارم 


در 9 همبستگی ضعيف و مثبت» در ۲ همبستكى متوسط و مثبت» در ۸ همبستكى قوى و مثبت» در 
7 همبستگی ضعيف و منفی» در ۶/ همبستكى متوسط و منفى و در 1٩‏ همبستگی قوی و منفى بين 
ميزبان فروندٌ ت و افت سطح آب زيرزمينى وجود دارد. مقادير همبستكى پیچیدگی ارتباط ميزان 


دهنده لزوم انجام بررسی‌های دقیق تر و ارائه مدل‌های جامع تر است. 


کلیدواژه‌ها: فرونشست. تداخل‌سنجی راداری. تحلیل سری زمانی» سطح آب زیرزمینی جاه‌های پیزومتری 


١‏ مقدمه 


افت تدريجى ارتفاع سطح زمين كه با جابه‌جایی‌های ناجيز افقى همراه است را فرونش ست كويند. به دليل 
تدريجى بودن آن و ايجاد خسارات غیرقابل‌جبران» اين يديده را زلزله خاموش می‌نامند. آب حاصل از بارندكى يا 
ذوب برفهاء در عمق خاک حركت می‌کند تا به لایه‌های غيرقابل نفوذ می‌رسد سپس تمام منافذ را ير می‌کند و 
منطقه اشباعى را به وجود می‌آورد كه به آن آبخوان (لايههاى آبدار يا آبخوانه) كفته می‌شود. برداشت آب زيرزمينى به 
علت افت سطح ايستابى و كاهش فشار سيال و افزايش فشار بين ذرات خاک منجر به ايجاد تراكم و يديده 
OPH PENET‏ 

oles‏ سیت هان مین نيك (GS) Be‏ ر شای ترازبالى Sa gos‏ براق بیش 
فرونشست است. هر دو اين روش‌ها محدود به ایستگاه‌های اندازه‌گیری هستند و نياز به حضور در منطقه انجام 


مشاهدات میدانی» صرف وقت و هزینه زياد دارند. از ديكر روش‌های ژئودتیکی يايش فرونة ت» تکنیک 


تداخلسنجى راداری است. به دليل تدريجى بودن يديده فرونشست برای بررسی آن بايد از روش‌های سرى زمانی 
استفاده کرد. ازجمله اين روش‌هاء تحليل سرى زمانی تداخلسنجی رادارى مبتنی بر پراکنشگرهای دائمى 
(PSINSAR)‏ و تحليل سرى زمانی تداخل سنجی راداری بر مبنای طول خط مبناى كوتاه (SBAS)‏ هستند. در اين 
تحقیق برای محاسبه فرونشست از روش PSINSAR‏ استفاده شده است» اين روش به دليل استفاده از نقاط پراکنشگر 
دائمی قادر است فاز ناشی از اتمسفر و عدم همبستگی زمانى و مکانی. خطاهاى باقیمانده توي وكرافى و نویز را کاهش 
دهد. در تحقیقات متعددی نتایج اين روش با مشاهدات GPS‏ و ترازیابی دقیق مقایسه شدهاند که نشان داده شده 


دقفت روش ۴ در حد میلی‌متر است (تامبورینی و همکاران. Ces‏ (روچی و همکاران. ۲( نتایج 


1 Global Positioning System 

2 Persistent Scatterer Interferometric Synthetic Aperture Radar 
3 small baseline subset 

4 Tamburini et al. 

5 Rucci et al. 


سال يازدهم تررسى ميزان. هميسكن :بين فرونشست و تغبيرات :+ ۱۷۵ 


اين روش همبسككى بالایی با مشاهدات GPS‏ و ترازيابى دقيق دارد (بيسواس و همكاران. SLE) OTA‏ و 
همکاران, "0٠6٠١‏ (خرمى و همكاران ۲۰۱۹)" تداخل‌سنجی راداری در دهه اخير به‌طور گسترده در زمينه مطالعات 
فرونش ست زمين موردتوجه پژوهشگران قرار گرفته است. در تحقيقى جابهجايىهاى پل Ka‏ كنك-زوهاى-ماكائو' 
با استفاده از روش AVL PSINSAR‏ تصوير ستتينل” از ٦‏ ژانویه ۲۰۱۸ تا ۲۷ نوامبر با استفاده از ۲۰۲۰ محاسبه شده 
است. نتایج اين تحقیق نشان داده است که حداکثر نرخ سالانه فرونشست در اين منطقه ۱۷ میلی‌متر است SL)‏ و 
همکاران ۲۰۲۱)'. در مقاله‌ای با استفاده از روش PS-InSAR‏ و داده‌های Cosmo-skymed‏ در بازه زمانی ۸ ژوئن 
۱ تا ۱۵ نوامبر ۲۰۱۷ تصاویر ستتینل در بازه زمانی ۱۸ دسامبر ۲۰۱۶ تا ۲۷ نوامبر ۲۰۱۸ و تصاویر 


ALOS ۴‏ در بازه زمانی ۱٩‏ ژانویه ۲۰۰۷ تا ۲۰ ژانویه ۲۰۱۱ جهت محاسبه فرونة ت منطقه دهلی ۲ 


استفاده شده است. در اين مقاله ميانكين جابه‌جایی‌های اين منطقه را با داده Cosmo-skymed‏ ۱۵ میلی‌متر در سال» 
حداکثر جابه‌جایی محاسبه شده از تصاویر ٠١ ALOS‏ تا ۱۸ میلی‌متر در سال و جابه‌جایی‌های حاصل از تصاویر 
ستتینل» ۲ تا ۱۳ میلی‌متر در سال آورده شده است (ملیک و همکاران2)۲۰۲۱ فرونشست شهر چانگچون جين را با 
استفاده از روش PSINSAR‏ با تصاویر سنتینل در بازه زمانی ۲۰۱۹ تا ۲۰۲۰ محاسبه شده و نتایج اين مقاله نشان داده 
است که حداکثر نرخ سالانه فرونشست اين منطقه. ۲۰/۶۵ میلی‌متر و ميانگین نرخ سالانه فرونشست ۲/۹۸ میلی‌متر 
استء همچنین در اين تحقیق همبستگی بين نتايج PSINSAR‏ و داده‌های ترازیابی درجه ۲ ۰/۷ به دست آمده که 
اين مقدار همبستگی به قابل‌اطمینان بودن نتایج PSINSAR‏ اشاره دارد (وانگ و همکاران, ۲۰۲۱)*. در تحقیقی دیگر 
نقشه فرونشست در جینگهای '' جين را با استفاده از روش PSINSAR‏ با تصاویر تراسار"" تهیه شده است» بیشترین 
مقدار فرونشست در سال ۲۰۰۹ در قسمت شرق منطقه مطالعاتی اتفاق افتاده كه اين مقدار بين ۶۰ تا ۵۰ میلی‌متر در 


rye aor‏ اسب (ليو و همکاران, ۰ فرونشست منطقه اولسان ۲ با روش PSInSAR‏ و SBAS‏ در بازه ۹ ت 


۱ بررسی کرده‌اند نتایج اين مقاله نشان داده است كه حداکثر نرخ فرونشست به دست آمده از روش ٤/٤ SBAS‏ 


1 Biswas et al. 

2 Chang et al. 

3 Khorrami et al. 
4 Hong Kong—Zhuhai—Macao 
5 Sentinel-1A 

6 Xiong et al. 

7 Dehli NCR 

8 Malik et al. 

9 Wang et al. 

10 Jinghai 

11 TerraSAR-X 
12 Liu et al. 

13 Ulsan 


۱۳ جغرافیا و مخاطرات محیطی شمارةً چهارم 


سانتی‌متر در سال و از روش PSINSAR‏ ۶ سانتیمتر است (جو و همکاران ۲۰۱۱). نرخ فرونشست در مکزیک در 
بازه زمانی ۲۰۰۳ تا ۲۰۰۷ با داده‌های ماهواره انویست" با استفاده از روش ۸۰۳۹105۸1 میلی‌متر در سال محاسبه 
کر ده‌اند. نتایج نشان داده که حداکثر فرونشست در اطراف جاه‌ها النت (بکسی و همکاران. ۹ تحقیقات متفاوتی 


جهت بررسی روابط بين فرونشست و عوامل تأثیرگذار بر OF‏ انجام شده است. در بعضی از اين تحقیقات تنها تأثیر عامل 


بررسی و بدون توجه به تاثيرات ساير عوامل» بين مقادیر فرو 


تغييرات سطح آب زيرزمينى بر يديده فرو و 
تغییرات سطح آب زيرزمينى رابطه خطی ارائه شده است. در تحقيقى نرخ فرونشست در بندر اینچتون * واقع در بوسان* 
o S‏ را با JERS-1°‏ با استفاده از روش PSINSAR‏ ۳۰ میلی‌متر در سال و در منطقه ایستگاه راه‌آهن ژوان" ۱۵ 
میلی‌متر در سال تخمین زده‌اند. درنهایت در اين مقاله دليل اصلی فرونشست در اين منطقه ساحلی را افزایش سطح 
آب دریا دانسته‌اند (كيم و همکاران, ۸6۲۰۰۷ در پژوهشی دیگر نقشه فرونشست دشت لیائوها" در جين را با تصاویر 
۴ با استفاده از روش 29125.41 در بازه ۲۰۰۷ تا ۲۰۱۱ تهیه شده است. طبق نقشه ارائه شده در 
اين مقاله حداکثر نرخ فرونشست در اين منطقه ۲۳۷ میلی‌متر در سال و كل مقدار فرونشست در ٤‏ سالء ۷۹۹ میلی‌متر 
است. همچنین در اين تحقیق نتیجه‌گیری شده است که بين فرونشست و مناطق استخراج نفت همبستگی وجود دارد 
و کشف شده که مقادیر بالای فرونشست TSI‏ در مناطقی است که چاه‌های نفت وجود دارند (سان و همکاران 


07 همچنین تحقیقی ارائه شده است که حداکثر نرخ فرونشست در شهر ye‏ ايتاليا را با دو مجموعه داده‌ی 
"ERS 1/2‏ و انویست با استفاده از روش PSINSAR‏ در بازه زمانی ۱۹۹۲ تا ۰۲۰۱۰ ۲۰ میلی‌متر در سال برای هر دو 
داده به دست lea yg!‏ مطالعاتی که روی بخش شرق اين منطقه انجام دادهاند نان داده است که بين قدمت 
ساختمان‌ها و نرخ فرونشست همبستگی وجود دارد به اين صورت که در ساختمان‌های با قدمت کمتر نرخ 
فرونشست کمتری تخمین زده شده است (سولاری و همکاران ۲۰۱۹) ". در تحقیقی دیگر فرونشست در دشت 
بئیجینگ *" جين در بازه زمانی ۲۰۱۱ تا ۲۰۱۸ با استفاده از روش ۳۹105۸ به دست آورده شده و سپس ارتباط بين 


RE‏ و تغییرات سطح آب زیرزمینی را در سال ۲۰۱6 و ۲۰۱۸ در چهار آبخوان بررسی کرده‌اند. حداکثر مقادیر 


1 Jo et al. 

2 Envisat SAR 

3 Békési et al. 

4 Incheon 

5 Busan 

6 Japanese Earth Resource Satellite-1 
7 Juan 

8 Kim et al. 

9 Liaohe 

10 Sun et al 

11 Pisa 

12 European Remote-Sensing Satellite 
13 Solari et al. 

14 Beijing 


سال يازدهم تررسی ميزان. همستگی ین فرونشست و تغبيزات ی ۱۷ 


کرولیشن ارائه شده در اين مقاله براى اين چهار آبخوان, ۰۵ eny‏ ۰/۸ و ۰/۵ انیت (جن و همکاران. A (Yeye‏ در 


مقاله‌ای فرونشست تایوان با جی‌پی‌اس محاسبه و سپس رابطه خطی بين فرونشست و تغییرات سطح آب زیرزمینی 
در بازه زمانى ۶ تا ۲۰۰۲ را به دست أوردهاندك در اين مقاله بر اساس نتايج» برداشت بی‌رویه از منابع آب زیرزمینی 
را تأثیرگذارترین عامل در رخداد فرونشسست اين منطقه دانسته‌اند (جن و همکاران ۲۰۱۰). فرونشست سمارتكى" را 
در بازه زمانی ۲۰۱۶ تا ۲۰۱۷ را با استفاده از روش PSINSAR‏ و تصاویر سنتینل محاسبه شده که حداکثر فرونشست 
در منطقه 0- سانتیمتر است» سپس همبستگی بين فرونشست و تغییرات سطح آب زیرزمینی را 4۷/۵۸ درصد به 


دست آورده‌اند (فیرداوس و همکاران 0018 " در پژوهشی رابطه بين فرونشست و تغییرات سطح آب زیرزمینی را 


را با استفاده از روش تداخل سنجی راداری محاسبه و سپس 
همبستگی بين نرخ فرونشست و نرخ افت سطح آب زیرزمینی را در سال ۲۰۰5 تا ۲۰۱۱ را محاسبه کرده‌اند (مقادیر 
همبستگی در اين بازه زمانی بين ۰/۱ تا 0/4 و رابطه خطى بين اين دو پارامتر را به دست آورده‌اند (زارعی و 
همکاران, ۱۳۹۹). بر خلاف مقالات ذکر شده که رابطه خطى و مستقیم بين فرونشست و تغییرات سطح آب زیرزمینی 
ارائه داده‌اند. نتایج اين تحقیق در مورد رابطه فرونشست و تغییرات سطح آب زیرزمینی در منطقه اصفهان نشان 
می‌دهد که در همه مناطق نمی‌توان رابطه خطى بين فرونة ت و سطح آب زیرزمینی به دست آورد و مقادیر 
همبستگی بين اين دو پارامتر با توجه به ویژگی‌های منطقه تغییر می‌کند. به دلیل اينكه فرونشست تحت تأثیر چندین 
عامل طبیعی و انسانی است نباید تنها اثر يك عامل بررسی شود. 


بررسی شده است. در اين تحقیق ابتدا فرو 


۲- منطقه موردمطالعه 

استان اصفهان در مرکز ایران واقع شده است. اين استان با مساحتی حدود ۱۰۵۹۳۷ کیلومترمربع بين ۳۰ درجه و 
Guay‏ ۶ د ىكل هه خرن بیان عط اراو ا gas‏ وا De) ypa‏ 
شرقی از نصف النهار كرينويج قرار دارد. بر اساس آخرین تقسیمات كشورىء اين استان دارای ۱۷ شهرستان» 7۰ شهر 


فرو 
باعث تشدید of‏ شده است. حداکثر درجه cule‏ در تابستان ۳۹ درجه gle‏ كراة است که تابستان‌هایی گرم و 
کاک ا زد وکر Bi Ba‏ ريات هرهم رسن وان یکی A‏ ارعان Toa réel‏ 
بناهای تاریخی متعدد دارای شهرت جهانى است. برای حفظ و مراقبت از بناهای تاریخی و به دلیل اهمیت و وسعت 
و جمعیت اين منطقه. بررسی و مطالعه‌ی فرونشست و جلوگیری از پیشروی آن در مناطق مسکونی حیاتی است. البته 


است که دو دلیل مهم تأثیرگذار عدم مدیریت صحیح منابع آب زیرزمینی وعدم جریان دائمی زاینده‌رود 


1 Chen et al. 
2 Semarang City, Indonesia 
3 Firdaus et al. 


\VA‏ جغرافیا و مخاطرات محیطی شمارةً چهارم 


در نزدیکی آثار تاریخی وجود دارد و در آینده آثار تاریخی در 


AS SX) اسن‎ 


فرسال‌های اتر مشاهداتی از وجود فقوو 


۳- داده‌های مورداستفاده 

داده‌های مورداستفاده در اين تحقیق شامل ۷۳ تصوير راداری سنتینل با فرمت SLC‏ در بازه زمانی ۲۰۱۶ تا ۲۰۲۰ 
در عبور بالاگذر" در مسیر" شماره ۲۸ ماهواره است اين تصاوير در باند C‏ معادل طول موج ۵/۵۶ سانتی‌متر و 
فرکانس ۵/4۰۵ KS‏ هرتزء با زاویه فرود ۳۹/۵۰۵۸ با دو پلاریزاسیون ۷۷,۷11 درمد 1۷۷ "و عرض تصویربرداری 
۰ كيلومتر با قدرت تفکیک مکانی ۳/۱۷۱۳/۹۵ متر اخذ شده‌اند. از مدل ارتفاعی رقومی ۹181 با قدرت تفکیک 


مکانی سه ثانيه ٩۰(‏ متر) در روند پردازش PSINSAR‏ برای محاسبه فرونشست استفاده شده است. 


حدول ۱- مشخصات تصاویر راداری مورداستفاده در پژوهش 


قدرت تفکیک در 
قدرت تفکیک در ۲ auge‏ بازه زمانی تعداد 
جهت ازیموت شماره سیر | مدتصویربرداری فرمت داده داده 
جهت رنج (متر) مدار تصاویر تصاویر 
(متر) 
UL ۱۳/۳۹۵ ۳/۳۷‏ گذر SLC IW YA‏ ۲۰۱۶-۶۰ ۷۳ ستتینل ۱ 


1 Ascending 
2 Track 
3 Interferometric Wide 


سال يازدهم بررسى ميزان همبستكى بين فرونشست و تغييرات ... \v4‏ 


51-300 ع‎ 51°45 0"E 


51° 15'0" 51° 30'0" 51°45'0°E 52°0°O°"E 


شکل ۱- محدوده موردمطالعه پژوهش 


برای حذف اثرات فصلی در همبستگی بين فرونشست و تغییرات سطح آب زیرزمینی در اين تحقیق تغییرات 
سطح آب زیرزمینی آبخوان‌ها و فرونشست در ماه شهریور سال‌های ۱۳۹۳ تا ۱۳۹۷ استفاده شده است. در شهریور 
ماه مقدار بارندگی منطقه موردمطالعه و برداشت از آبخوان حداقل است. بنابراین اين ماه بهترین ماه برای بررسی 
سطح آب زیرزمینی است. از نقشه زمین‌شناسی با مقیاس ۰ منطقه اصفهان برای بررسی خحصوصیات جنس 
خاک منطقه و تأثير آن بر فرونشست استفاده شد. با توجه به اين نقشه جنس خاک در منطقه موردمطالعه در اين 
پژوهش آبرفت» شیل» تراس و سنگآهک است. تمام چاه‌های پیزومتری مورداستفاده در اين تحقیق در مناطق با 


جس آبرفت و تراس واقع شده‌اند. 


۱۸۰ جغرافیا و مخاطرات محیطی شمارةً چهارم 


4- تداخل سنجی راداری 

تداخل‌سنجی راداری تکنیکی است که از تداخل امواج الکترومغناطیس جهت رسیدن به اطلاعات مفید مانند تهیه 
مدل رقومی ارتفاعی» میزان جابه‌جایی زمين و ... استفاده می‌کند. اين تکنیک از اختلاف فاز زوج تصوير راداری 
استفاده می كند, بنابراین پارامتر فاز امواج نقش اساسی را در اين تکنیک دارند. روش تداخل‌سنجی تفاضلی مرسوم و 
روش‌های تحلیل سری زمانی مانند طول خط مبنای کوتاه" و تحلیل سری زمانی تداخل سنجی راداری مبتنی بر 
پراکنشگرهای دائمی از روش‌های تداحل‌سنجی تفاضلی است. با وجود اينكه روش تداخل‌سنجی تفاضلی مرسوم" 
جابه‌جایی‌های زمین را با دقت خوبی اندازه‌گیری می‌کند. نمی‌تواند نقاطی را که بازپراکنش آن‌ها با زمان تغییر می‌کند 
به دليل عدم همبستگی‌های زمانی موردبررسی قرار دهد. مناطق پوشیده از Gy‏ زمین‌های با پوشش ALS‏ و 
جنگل‌ها ازجمله مناطق با مشکل عدم همبستگی هستند. با توجه به اينكه در روش‌های تحلیل سری زمانی از 
مجموعه‌ای از تداخل‌نماها در يك بازه زمانی استفاده می‌شود عدم همبستگی‌های مکانی و زمانی و اثرات اتمسفر 
کاهش می‌یابد بنابراین اين دسته از روش‌ها از روش تداخل سنجی راداری تفاضلی مرسوم دقیق‌تر هستند. 

۱-۶-تحلیل سری زمانی تداخل سنجی راداری مبتنی بر پراکنشگرهای دائم PSINSAR))‏ 

ازجمله محدودیت‌های روش تداخل‌سنجی راداری مرسوم. عدم همدوسی تصاویر مانع از تولید تداخل‌نماهای 
تفاضلی می‌شود. حطای توپوگرافی باقیمانده در تداخل‌نما (بزرگی اين حطا به کیفیت DEM‏ و خط مبنای قائم وابسته 
است) می‌باشد که اين محدودیت‌ها سبب توسعه روش PSINSAR‏ شده است. روش PSINSAR‏ استاندارد را اولین 
بار فرتی و همکاران (۲۰۰۱) مطرح کردند اين روش با هدف غلبه بر محدودیت‌های تکنیک تداخل‌سنجی معمول 
pte)‏ همبستگی‌های مکانی و زمانی و اثر اتمسفر) با استفاده از نقاط همدوس در بازه زمانی آنالیز به‌دقت پایش 
میلی‌متری فرونشست می‌رسد. درصورتی که تراکم نقاط پراکنشگر کم باشد. تخمین فاز اتمسفر و حذف آن امکان‌پذیر 
نیست (فرتی و همکاران. 6۲۰۰۱" روش ۳۹108۸1 به‌عنوان یک ابزار قوی برای تخمین جابه‌جایی‌های کوچک در 
مناطق شهری بزرگمقیاس در طول بازه‌ی زمانی‌های بلندمدت يا کوتاه‌مدت شناخته شده است. پراکنشگرهای دائمی 
نقاطی است که در طول زمان فاز آن‌ها ثابت می‌ماند» با استفاده از اين پراکنشگرهای دائمی فاز ناشی از اتمسفر و 
عدم همبستگی زمانى و مکانی. خطاهاى باقيمانده تويوكرافى و نويز حذف می‌شود. اشكال این روش فقدان پیوستگی 
بین داده‌هاست» داده‌ها مجموعه‌ای از نقاط است که پراکندگی آن‌ها به شكل و کاربری منطقه بستگی دارد. اين نقاط 


1 Differential Interferometric Synthetic Aperture radar (D-InSAR) 
2 Small Baseline Subset (SBAS) 

3 Differential Interferometric Synthetic Aperture radar (D-InSAR) 
4 Ferretti et al. 


سال يازدهم تررسى ميزان. هميسكن :بين فرونشست و تخیرات :+ ۱۸۱ 


PAR‏ بهترين عملكرد را در مناطق شهرى و مناطق با جنس سنگ دارد (اوستير و همکاران, ۲۰۰۷). داده‌های 
موردنيازء 2+1 تصوير راداری منطقه با فاصله زمانىهاى مناسب است كه یک تصوير به‌عنوان تصوير يايه (master)‏ و 
2 تصوير به‌عنوان تصاوير پیرو (Slave)‏ در نظر گرفته می‌شوند. تصوير ah‏ بر اين اساس انتخاب مىشود كه عدم 
همبستگی مكانى و زمانى تصاوير مينيمم شود به‌عبارت‌دیگر به ازاى تداخل نماهاى تشكيل شده برمبناى آن» مقدار 
همدوسى تجمعى Ol‏ بيشينه شود. بعد از انتخاب نقاط يراكنشكر دائمى در مجموعه تداخل‌نماها با یکی از سه روش 
پراکندگی دامنه بر مبناى شاخحص پراکندگی دامنه. روش مبتنی بر پایداری فاز و روش مبتنى بر همبسستگی, بين اين 
نقاط شبکه‌ای با استفاده از روش مثلث‌بندی دلونی " تشکیل می‌شود. اين شبکه اساس تخمین فاز ناشی از اتمسفر و 
افزايش نقاط PS‏ است. به علت مجاور بودن نقاط در شبکه و همبسته بودن اتمسفر در مکان» اختلاف فاز بين دو 
پیکسل همسایه فاز اتمسفر حذف خواهد شد (هانسن" ۲۰۰۱). در مثلث‌بندی بين نقاط اگر طول اتصالی از 
همبستگی اتمسفر بیشتر باشد آن اتصال حذف می‌شود تا افزونگی اتصالات کاهش LL‏ و دقت تخمین جابه‌جابی 


زمين افزایش یابد. 


۵ نتایج تحقيق 

به‌منظور بررسی جابه‌جایی‌های زمین در اصفهان و چند منطقه اطراف آن, نقشه جابه‌جایی با ۷۳ Sentinel- „ pai‏ 
۸ از روش PSINSAR‏ تهیه شد. حداکثر مقدار تجمعی فرونشست منطقه اصفهان در بازه ۲۰۱۶ تا ۲۰۲۰ تقريباً به 
۰ سانتی‌متر رسیده است؛ با توجه به نقشه زمین‌شناسی و نقشه فرونشست تجمعی منطقه. مشخص است که 
فرونشست در مناطق با جنس آبرفت و تراس اتفاق افتاده است و در مناطق با جنس شيل و سنگآهک فرونشست 
مشاهده نشده است (شکل AY‏ 


1 Oštir et al. 
2 Delaunay triangulation 
3 Hanssen 


۱۸۲ جغرافيا و مخاطرات محيطى شمارة چهارم 
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شکل ۲- مقادیر تحمعی فرونشست منطقه موردمطالعه در بازه ۶ تا ۲۰۲۰ 


سری زمانی فرونشست و تغیبرات سطح آب زیرزمینی در سال ۱۳۹۷ در چاه جاده پرزان در دشت مهيار جنوبی - 
دشت آسمان و كروج و شرودان دشت فلاورجان رسم شده است (شکل ۳ و شکل CE‏ همبستگی اين دو پارامتر در 
چاه جاده پرزان در طول یک سال ۰/۳۱- و در كروج و شرودان ۶ است که اين نشان می‌دهد در اين بازه یک ساله 
اين دو پارامتر رفتار مستقلی دارند. به ازاى افت سطح آب زیرزمینی در یک چاه ممکن است فرونشست با تأخير 
زمانی اتفاق بیفتد و دقیقاً همزمان با افت آب. فرونشست مشاهده نشود كه البته در اين مسئله جنس خاک منطقه در 
رخداد فرونشست تأثیرگذار است. نوسانات سطح آب در ماه‌های مختلف به دلیل برداشت‌ها و بارندگی نيز در اين 
همبستگی تأثیرگذار است بنابراین بازه زمانی يك سال برای بررسی رابطه فرونشست و تغییرات سطح آب زیرزمینی 


سال يازدهم تررسی ميزان همستگی :بين pind‏ تغيرات:: ۱۸۳ 


Land subsidence vs. groundwater expoitation 
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Land subsidence vs. groundwater expoitation 
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شكل -٤‏ سرى زمانی تغیبرات فرونشست و سطح آب زيرزمينى جاه كروج و شرودان در سال ۱۳۹۷ 


برای بررسى ميزان تغييرات سطح آب زيرزمينى چاه‌های پیزومتری نمودار تغييرات سطح آب زيرزمينى چند 
نمونه جاه در دشتهاى منطقه موردمطالعه رسم شد. نمودار تغييرات سطح آب زيرزمينى جاه ياقوت آباد جديد و 
جاه جاده جعفر آباد در محدوده اصنفهان-بر خوار و جاه لورک محمدیه و اراضی سده لنجان در دشت لنجانات در 


شكله و شكل 5 آورده شده است» نمودارها بيانكر روند كاهشى سطح آب زيرزمينى در اين چاه‌ها است. مقدار 


۱۸۶ جغرافیا و مخاطرات محیطی شمارة چهارم 


کاهش سطح آب زیرزمینی چاه جاده جعفرآباد ۲۰/۲۹- متر و در یاقوت abl‏ جدید ۳۱/۱۱- مت در لورک محمدیه 


7- متر و در اراضی سده لنجان ۳/۵۵- متر در بازه زمانی ۱۳۸۰ تا ۱۳۹۸ است. 
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شكل ۵- تغييرات سطح آب زيرزمينى جاه جاده جعفرآباد و ياقوت آباد جديد 
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شکل -٦‏ تغییرات سطح آب زیرزمینی لورک محمدیه و اراضی سده لنجان در دشت لنجانات 


برای مقایسه روند تغيير سطح آب زیرزمینی و فرونشست برای نمونه در محدوده اصفهان-برخوار چاه يونارت و 
چاه سه راهی سهل‌آباد. رفتار افت سطح آب زیرزمینی و رفتار فرونشست در بازه ۱۳۹۳ تا ۱۳۹۷ در شكل ۷ و شكل 


A‏ آورده شده است. نمودارها نشان می‌دهند كه در بازه ۱۳۹۳ تا ۱۳۹۷ در يونارت به ازای ۶/۸٩‏ متر افت سطح آب 


سال يازدهم بررسی میزان همبستگی بين فرونشست و تغییرات ... ۸۵ 
زیرزمش ۱١‏ سای قرو شنب و در سه راهی abl few‏ به ازای ۱/۱۳ متر افت سطح آب زیرزمینی» ۱۳ gaila‏ 
فروندٌ ت داشته است. مقدار 2 iua‏ فروندٌ نشست و افت سطح آب زیرزمینی در یونارت ۰/۹۹ و در سه راهى 


سهل‌آباد ۰/۹۵ است. 


Land subsidence vs. groundwater exploitation 
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شکل ۷- فرونشست و تغییرات سطح آب زیرزمینی در منطقه اصفهان برخوار (یونارت) 


Land subsidence vs. groundwater exploitation 


شکل ۸- فرونشست و تغییرات سطح آب زیرزمینی در منطقه اصفهان برخوار (سه راهی سهل آباد) 


در محدوده اصفهان کوهپایه سگزی (جوزدان) و محدوده فلاورجان (کروچ و شرودان) رفتار تغییرات سطح آب 
زیرزمینی و رفتار فرونٌ a‏ در بازه ۱۳۹۲ تا ۱۳۹۷ در شكل 9و شكل ٠‏ آورده شده استهء نمودارها نشان 


۱1۸٦‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


می‌دهند که در بازه ۱۳۹۳ تا ۱۳۹۷ در جوزدان به ازای ۵/۳۹ متر افت سطح آب زيرزمينى» ۸ سانتیمتر فرونشست» در 


كروج و شرودان به ازای ۱0/۲۵ متر افت آب زیرزمینی» ۷ سانتیمتر فرونشست داشته است. مقدار همبستگی 


Land subsidence vs. groundwater exploitation 
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شکل -٠١‏ فرونشست و تغییرات سطح آب زیرزمینی در كروج و شرودان 


در محدوده مهيار جنوبی-دشت آاسمان (شمال آجر تاکی) و محدوده ; نج فاباد (غرب درجه) نمودار خطی 
و ست و تغییرات سطح آب زیرزمینی در بازه ۱۳۹۳ تا ۱۳۹۷ رسم شده است (شکل ۱۱ و شکل VY‏ در شمال 


سال يازدهم تررسی ميزان. همستکی :بين فر pote iy‏ تغبيرات :: \AV‏ 
آجر تاكى به ازای ۱/۱۳ متر افت سطح آب زیرزمینی» ۲ سانتيمتر فرونشست و در غرب درچه به ازاى T‏ سانتيمتر 


افت سطح آب زیرزمینی» ۳ سانتیمتر فرونشست داشته است. مقدار همبستگی فرونشست و افت سطح آب زیرزمینی 
در منطقه شمال آجر تاکی ۰/۳۹ و در غرب درچه ۰/۰۸ است. 


Land subsidence vs. groundwater exploitation 


p = 5 


شکل ۱۱- فرونشست و تغييرات سطح آب زیرزمینی در منطقه مهيار جنوبی-دشت آسمان (شمال آجر تاکی) 


Land subsidence vs. groundwater exploitation 


p= 0.08 


شکل ۱۲- فرونشست و تغیبرات سطح آب زیرزمینی در منطقه نجف آباد (غرب درچه) 


تمام پیزومترها در منطقه موردمطالعه در اين تحقیق در مناطق با خاک آبرفت و تراس هستند به همین دلیل افت 


سطح آب زیرزمینی باعث می‌شود. آب از بين ذرات خاک برداشته شده و تنش مؤثر وارد بر ذرات خاک افزایش يافته 


مرا جغرافیا و مخاطرات محیطی شمارة چهارم 


و خاک متراكم مى شود درنتيجه فرونشست در مناطق با خاک آبرفتی و تراس مشاهده شده است ولى در مناطق با 
جنس شيل و سنگآهک به دليل تحكيم و تراکم‌پذیری پایین آن‌ها فرونشست در اين مناطق مشاهده نشده است. 
مقادير همبستكى فرونشست و تغييرات سطح آب زيرزمينى در بازه زمانى ۱۳۹۳ تا ۱۳۹۷ درشكل ۱۳ شكل ۱۳ 
آورده شده است. طبق اين نقشه ۱۲/ همبستگی‌های بين ۰/۲۵- تا ۰/۲۵ كه در اين بازه همبستگی بين دو يارامتر 
وجود ندارد بنابراين در اين مناطق علاوه بر افت سطح آب زيرزمينى عوامل ديكرى مثل جنس خاک شیب جهت 
شيب گسل و ... در رخداد فرونشست تأثیرگذار هستند. 19 همبستكى ضعيف و Cate‏ (همبستگی‌های بين ۰/۲۵ تا 
۰/۵ ۲ همبستكى متوسط و مثبت (همبستكى بين ۵ تا 0 و 1۸/ همبستگی قوى و مثبت (همبستكى بين 
۵ تا ۱) / همبستگی ضعیف و منفی (۰/۲۵۱-۰/۵- 4/ همبستگی متوسط و منفی (۰/۷۵- تا ۰/۵ و JA‏ 
همبستگی قوی و منفی (۱- تا ۰/۷۵-) بین فرونشست و تغییرات سطح آب زیرزمینی را نشان می‌دهند. محدوده‌ها با 
همبستگی قابل‌توجه مثبت (۰/۵ تا )١‏ به اين معنی است که در اين مناطق به ازای افت سطح آب زیرزمینی فرونشست 
در منطقه اتفاق افتاده است و جنس خاک منطقه هم آبرفت و رس بوده است که باعث شده است با تخلیه آب از 
منافذ خاک خاک متراکم شده و فرونشست اتفاق بیفتد. درصورتی که خاک منطقه از نوع كانى رس باشد بنابراین رس 
فعال است در صورت متراکم شدن غیرقابل نفوذ می‌شوند و آبخوان قابلیت احيا نخواهد داشت ولی درصورتی‌که از 
نظر دانهبندى» رس گفته شود ممکن است خاک ریزدانه باشد ولی از نوع کانی رس نباشد و غیرفعال باشد در اين 


صورت ممکن است قابلیت نفوذپذیری داشته باشد و دوباره آب در آن نفوذ کند. 
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شكل ۱۳- ضريب همبستگی بين فرونشست و تغييرات سطح آب زيرزمينى در بازه ۱۳۹۳ تا ۱۳۹۷ 


سال يازدهم بررسى ميزان همبستكى بين فرونشست و تغييرات ... ۱۸۹ 


1- نتيجه گیری و پيشنهادها 

در اين تحقیق تغییرات سطح آب زیرزمینی جاه‌های پیزومتری دشت‌های اصفهان-برخوار» لنجانات» مهيار جنوبى» 
فلاورجان. نجف آباد و كوهيايه سگزی در بازه زمانی ۰ تا ۱۳۹۸ و همچنین همبستگی بين فرونشست و تغییرات 
سطح آب زیرزمینی در بازه زمانی ۱۳۹۳ تا ۱۳۹۷ در منطقه اصفهان آورده شد. مقادیر همبستگی بين اين دو پارامتر نشان 
می‌دهد که با توجه به اينكه عوامل مختلفی مثل جنس خاک منطقه. شیب جهت شیب ارتفاع و تغییرات سطح 


آب‌های زیرزمینی و ... در رخداد پدیده فرونشست دخالت دارند» بررسی تنها یک عامل روش صحبحی برای پی 


بردن به علل به وجود آمدن فرو ت نیست. ارائه روابط بين تغییرات سطح آب زیرزمینی و فرونشست ممکن 
است برای یک منطقه حاص قابل قبول باشد Jy‏ در حالت IS‏ بايد خصوصیات ویژه هر منطقه مثل جنس خاک 
كسلء شیب و.. هم در كنار تغييرات سطح أب زیرزمینی در نظر گرفته شود برای مثال در منطقه‌ای که جنس خاک؛ 
خاک ریزدانه باشد. برداشت بی‌رویه از منابع آب زیرزمینی موجب فرونشست می‌شود و ممکن است همبستگی بين 
اين دو پارامتر UL‏ و رابطه خطى به دست آید ولی در منطقه‌ای با جنس خاک درشت دانه» ممکن است در صورت 
افت زياد سطح آب زیرزمینی فرونشست کمتری مشاهده شود. رابطه خطی بين تغیبرات سطح آب زیرزمینی و 


فرونشست زمین فقط در حالت تثوری امکان‌پذیر است ولی در واقعیت رابطه فرونشست با چندین عامل مؤثر بر آن 


پیچیده‌تر از رابطه خطی است و عوامل مختلف مؤثر بر فرونشست بايد همزمان بررسی شوند. تأخير زمانی بين افت 
سطح آب زیرزمینی و رخداد فرونشست هم در همبستگی بين اين دو پارامتر مؤثر است. بنابراین رابطه فرونشست و 
تغییرات سطح آب زیرزمینی بايد در بازه زمانی طولانی‌مدت بررسی شود. Gb‏ نتایج حاصل از اين تحقیق, در مناطقی 
که همبستگی بين فرونشست و تغییرات سطح آب زیرزمینی ضعیف است يا همبستگی وجود ندارد عوامل تأثي ركذار 
دیگر بايد موردبررسی قرار گیرند تا با احتمال بیشتری بتوان دلیل فرونشست در اين مناطق را شناسایی کرد. همان‌طور 
که در نتایج هم آورده شده است در بعضی از اين مناطق با وجود افزايش سطح آب همچنان فرونشست زمین ادامه 
داشته است که اين می‌تواند نشان‌دهنده‌ی عدم برگشت‌پذیر بودن منطقه نسبت به فرونشست باشد به‌عبارت‌دیگر 
حتی با احیای مجدد منابع آب زیرزمینی در اين مناطق نمی‌توان شاهد بهبود وضعیت فرونشست بود. در اين موارد 
yok satel‏ عياف و کانی‌شناسی انجام شوه به‌طور مقال فرصورتی که خاک مطقه از نظر كات dlos a‏ 
اين رس فعال است بنابراین در صورت برداشت بی‌رویه از آب زیرزمینی خاک متراکم شده و قابلیت جذب آب را 
نخواهد داشت و فرونشست رخ می‌دهد. درصورتی که از نظر دانه‌بندی, منطقه از رس باشد. اين رس غیرفعال است 
و نسبت به خاک رس فعالء قابلیت جذب آب بیشتری دارد. در اين تحقیق نشان داده شد که بررسی رابطه 
فرونشست و تنها يك عامل تأثیرگذار بر آن» صحیح نیست و در یافتن علت رخداد فرونشست کمکی نمی‌کند. زیرا 
پدیده فرونشست همزمان ممکن است تحت تأثير چند عامل در منطقه موردنظر به وجود آید. پيشنهاد می‌شود برای 


۱۹۰ جغرافیا و مخاطرات محيطى شمارة چهارم 


خاک منطقه استفاده و همچنین تأثیر عوامل دیگر مثل گسل و شیب ارتفاع منطقه و ... همزمان بررسی شوند. 
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